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ABSTRACT 

> 

£ — ' The connection between some features of the mctallicity gradient in the Galactic 

C*~) , disk, best revealed by Open Clusters and Cepheids, and the spiral structure, is ex- 

plored. The step-like abrupt decrease in metallicity at 8.5 kpc (with Ro— 7.5 kpc, or 
C*~) ' at 9.5 kpc if Ro = 8.5 kpc is adopted) is well explained by the corotation ring-shaped 

. gap in the density of gas, which isolates the internal and external regions of the disk 

C^) ' one from the other. This solves a long standing problem of understanding the different 

chemical characteristics of the inner and outer parts of the disk. The time required to 
build up the metallicity difference between the two sides of the step is a measure of 
the minimal life-time of the present grand-design spiral pattern structure, of the order 
of 3 Gyr. The plateaux observed on each side of the step are interpreted in terms of 
the large scale radial motion of the stars and of the gas flow induced by the spiral 
structure. The star-formation rate revealed by the density of open clusters is maxi- 
mum in the Galactic radial range from 6 to 12 kpc (with an exception of a narrow gap 
at corotation), coinciding with the region where the 4-arms mode is allowed to exist. 
We argue that most of the old open clusters situated at large galactocentric radii were 
born in this inner region where conditions more favorable to star-formation are found. 
The ratio of a-elements to Fe of the sample of Cepheids does not vary appreciably 
with the Galactic radius, which reveals an homogeneous history of star formation. Dif- 
ferent arguments are given showing that usual approximations of chemical evolution 
models, which assume fast mixing of metallicity in the azimuthal direction and ignore 
the existence of the spiral arms, are a poor ones. 

Key words: Galaxy: Abundances, Galaxy: Evolution, Galaxy: Structure, Stellar 
Dynamics 



1 INTRODUCTION 

The variation of the abundance of chemical elements as a 
function of radius in the Galactic disk provides constraints 
to star formation history in the disk and to basic assump- 
* E-mail: jacques@astro.iag.usp.br tions made in chemical evolution models. Several reviews 
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have been published recently, focusing in different chemical 
elements like Iron or a-elements, and in different tracers like 
Open Clusters, Cepheids, H1I regions, etc., see e.g. the se- 
ries of p apers by Andri e vsky a nd collaborator s from 2002 
to 2005. iRudolph et all l|2006l ), iPedicelli et "all (|2009l) , and 
iMagrini et al l (|2009h . among others. In the past, the metal- 
licity gradients were most often fitted by straight lines in 
logarithm of abundance versus Galactic radius plots, with 
slopes of about -0.07 dex/kpc. This "first order" approxi- 
mation did not pose severe constraints to chemical evolu- 
tion models, since many important parameters of the mod- 
els, like gas density and infall rate of matter are supposed 
to be smoothly decrea sing functions of the Galactic radius 
l|Chiappini et al.lll997l ). However, evidences have accumu- 
lated that a single slope approximation (or two slopes, in 
the last reference) is not satisfactory, since it hides more 
complex slope changes that are real and should be better un- 
derstood. For instance, the series of papers reporting obser- 
vations of Cepheids by Andrievsky and collaborators showed 
that the slope of the abundance gradient is larger than av- 
erage around 5-6 kpc from the center, followed by a plateau 
from about 6.5 kpc to 10 kpc, and then by a new decrease 
at about 10 kpc (here we quote the Galactic radii as given 
in the original papers). As another example. iTwarog et al.l 
(1 19971 ) (hereafter TAAT) analyzing a sample of Open Clus- 
ters, have found that the Galactic metallicity gradient can 
be described in terms of two zones with flat slopes: an inner 
disk with radii ranging from 6 to 10 kpc and an outer disk 
for radii larger than 10 kpc. These two zones are separated 
by a step-like disco n tinuit y of about 0.3 dex i n [Fe/ Hj. Fi- 
nally, ICarnev et ail (|2005l ) and ICarraro et"ai] (|2007f ) found 
that the metallicity gradient is totally flat in the outer disk. 
All these results are in conflict with the old view of an unique 
and constant slope, however, they are not inconsistent with 
each other, as it will be shown in this work. 

These changes of slope and step-like behavior of the 
gradients merit great attention since they provide a much 
deeper understanding of the disk evolution and impose real 
constraints to the existing models. In particular, we believe 
that the effect of the spiral arms cannot be ignored in chem- 
ical evolution models, as it has most often been done. The 
spiral arms are the star-forming machines of the disk, and 
are the places where the youngest stars (O, B stars and 
HII regions) are found. The massive stars are those which 
have the most important role in chemical evolution. Fur- 
thermore, when one integrates the orbits of Open Clusters 
in the Galactic potential for a time equal t o their age, one 
finds that they were born in the spiral arms l|Dias fc Lepind 
2005). This is a confirmation that almost all stars are born in 
the arms. The spiral arms obey a number of rules dictated 
by the resonances between the epicycle frequency k (the 
frequency of oscillation of a star around its non-perturbed 
circular orbit) and the frequency of rotation u around the 
Galactic center, as observed in the frame of reference of the 
spiral pattern uj p . Among the major resonances are the in- 
ner and the outer Lindblad resonances (ILR and OLR, or 
2:1 resonances), situated where uu — uj p — ±«/2. According 
to the clas s ical t heory of galacti c spira l waves proposed by 
iLin fc Shul || 19641 ) and lLin et ail (|l969l ). the spiral arms are 
restricted to the interval between these two resonances, so 
that an important change in the star-formation rate could 
be expected at the corresponding radii. For other theories 



in which the spiral arms are explained by the crowding of 
stellar orbits, the ILR and OLR have the same nature and 
are equally important (|Kamaislll973l ). The 4:1 resonances, 
which appear where ui — cu p = ±k/4, are the equivalent of 
the ILR and OLR for 4 arms; in principle a 4-arms structure 
could only exist between these two resonances. A stellar or- 
bit at a 4:1 resonance presents 4 maxima of the radius during 
one revolution around the center, which gives it the aspect 
of a square orbit with rounded corners. 

Another major resonance is corotation, situated where 
the rotation speed of the material of the disk coincides with 
that of the spiral pattern (w = u) p ). This maximizes the 
effect of the potential perturbation of the spiral arms on the 
stars and gas clouds in this region. The spiral arms are able 
to produce systematic radial transfer of interstellar gas with 
opposite directions on the two sides of corotation, and to 
create a ring with a void of gas at this resonance, as it will 
be discussed later in this paper. 

A matter of concern is that it has been pointed out 
by several authors who analyzed the metallicity gradients in 
the Galactic disk that deviations of individual measurements 
from the average slope are often found to be larger than in- 
dividual errors on t he measurements (e.g. lYong et al.ll2005l . 
ICarraro et ai1l2007l) . This challenges the usually accepted 
hypothesis that the metallicity must have a single value 
at any given radius (or, equivalently, that the azimuthal 
diffusion is very fast). It is important to analyze the non 
axis-symmetric variations of metallicity, and if they are con- 
firmed, to verify if they are related to the spiral structure, 
or alternatively, if random local variations in metallicity do 
exist. 

In the present paper we analyze the connection of vari- 
ations in the abundances of elements with resonances of the 
spiral structure like the 4:1 resonance, corotation and OLR. 
The radii of these resonances, at least for the main spiral 
structure, is no more uncertain as they were in the past; 
they are now known with accura cy of the order of 0.2 kpc, as 
discussed bv lLepine etafl (|2010l ). hereafter L+5. We adopt, 
for corotation, R c — 8.4 kpc, combined with a solar orbit 
with Ro = 7.5 kpc. The reasons for the choice the "short 
scale" for Ro are discussed in Section 3. The value of R c re- 
sults from the nice agreement between parameters directly 
determined, such as the rotation speed of the arm pattern, 
measured by integration of the orbits of Open Clusters, the 
radius of the ring-shaped gap in the Galactic disk (associ- 
ated with corotation), and the observation of the position of 
the 4:1 resonance. One can find in the literature a number of 
determinations of the corotation radius that are in conflict 
with these measurements. However, we must be conscious 
that numerical N-body simulations, hydrodynamical simu- 
lations and chemical evolution models depend on a series of 
assumptions of unknown parameters and on simplified hy- 
potheses, so that their results cannot be regarded as being 
as precise as th ose of direct measurem ents. For instance, in a 
recent work bv lAcharova et all (|2010l ). a chemical evolution 
model gave a best fit to the metallicity plateau of Cepheids 
with corotation at 7 kpc (which would be about 6.6 kpc in 
the shorter scale adopted here). However, it was left clear in 
the mentioned paper that this result could not be regarded 
as a determination of the corotation radius. 

We are aware that dynamical measurement of the pat- 
tern speed of the spiral structure, based on the location 
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of the Hyades and Sirius moving groups in velocity space 
l|Famaev et al.|[2 005, 2007 ), placed the corotation at a radiu s 
R=12 kpc ijQuillen fc Minchevll2005l . IPompeia et"ai1l2011n . 
This could perhaps be reconciled with the present findings 
in the case of the existence of multiple spiral patterns with 
different pattern speeds; while the main grand-design spiral 
pattern has its corotation at 8.4 kpc, and is responsible for 
the plateau in the metallicity gradient as well as for the ring- 
shaped gap in the density of gas, it could coexist with an 
outer m=2 pattern whose 4:1 inner resonance smoothly con- 
nects with the corotation of the main pattern. This i s wha t 
happens in the N-body simulations of iQuillen et al.l (2010). 
While the corotation of the main spiral pattern would not 
have a large effect on the local velocity distribution of stars 
in the solar neighbourhood, the 4:1 inner resonance of that 
slower outer pattern would be responsible for the observed 
substructures. The corotation of this outer spiral pattern 
would then be located at about R = 12 kpc, near the OLR 
of the main structure, and could also contribute to a gap 
which seems to exist in the distribution of clusters at this 
radius, and will be discussed in this paper. We restate, how- 
ever, that although we wish to keep open the possibility of a 
second pattern speed, all the analysis in this paper is based 
on what we call the main structure, with corotation at 8.4 
kpc. 

The accuracy now reached by direct measurements of 
the radii of resonances of the main structure allows us to 
perform detailed comparisons between the features of abun- 
dance variations and the Galactic spiral structure. This is 
the main focus of the present paper, which is organized as 
follows. In Section 2 we analyze the expected ranges of radii 
covered by orbits of stars, to estimate how far from their 
birthplace young tracers can be found. This is a theoretical 
background needed for a better understanding of the abun- 
dance plateaux that are discussed later in the paper. In the 
subsequent two sections we present the data of the abun- 
dance tracers that have the most precise distances: Open 
Clusters (Section 3) and Cepheids (Section 4). In Section 
5 we examine a departure from axis-symmetry in chemical 
abundance that we call azimuthal gradient, and we comment 
an extra-galactic example of non-axissymmetric gradients 
which looks like the "overlapping" gradients that we believe 
to exist in the Galaxy. In Section 7, from the observation of 
the step-like decrease of metallicity at corotation, we derive 
a lower limit to the life-time of the spiral structure. 



2 RANGE OF RADII SWEPT BY STARS 

We review here the basic rules that regulates the range of 
radii swept by the orbits of stars or Open Clusters in the 
Galaxy. 

A typical orbit is shown in Figure 1. This orbit was ob- 
tained by numerical integration for a star with initial veloc- 
ity 40 km/s in excess with respect to circular rotation at the 
initial radius. If we are interested on the radial motion only, 
there is a classical solution which consists in introducing an 
effective potential governing the radial motion of a particle 
of a given angular momentum £ in a, central force field. The 
effective potential is the sum of the original central potential 
plus a rotational energy given by U — £ 2 /2mr 2 , where m is 
the mass of the particle and r its distance from the center. 



20 



10 : 



r 



-10 




-20 



-20 



-10 





X - kpc 



10 



20 



Figure 1. Typical orbit of a star in the Galactic potential, with 
minimum radius of about 7 kpc and maximum of about 11 kpc. 



Figure 2 shows two examples of effective potential curves, 
with different angular momentum of the star. The Galactic 
potential is obtained by integration of the rotation curve of 
the Galaxy (the gravitational force is equal to mV(r) 2 /r, 
where V(r) is the rotation curve). We arbitrarily set the po- 
tential equal to zero at a large distance (we used 60 kpc). 
Note that adding an arbitrary constant to the potential en- 
ergy curves does not change the dynamics of the stars, since 
we are not investigating the escape probability and we re- 
strict our study to orbits with small perturbation energy. 
In the figure, a dot represents the radius of minimum en- 
ergy (circular orbit), and the dashed lines show the range of 
radi us swept by the sta r for a not too large perturbation. 

iLepine et al.l (|200Sl ) presented an histogram of the per- 
turbation velocities of the Open Clusters (the residual veloc- 
ity after subtracting the velocity of the rotation curve). The 
histogram presents a peak at about 15 km/s, and extends 
to about 50 km/s. The ranges of radii illustrated in Figure 
2 were computed with a perturbation velocity of 40 km/s, 
so that they are close to the maximum observed range. The 
most important result is that the stars travel between their 
normal maximum and minimum Galactic radius in a time 
shorter than one revolution around the Galactic center. An- 
other mechanism capable of producing changes in the Galac- 
tic radius of stars is the scattering by the corotation reso- 
nance . Thi s was investigated in detail bv lSellwood fc BinnevI 
(2002) and lLepine et al] (2003), showing that the scattering 
is able to impart larger radial amplitudes, but it requires a 
longer time to be effective, compared to the normal epicycle 
oscillations described previously. 



3 GRADIENTS BASED ON OPEN CLUSTERS 

Our analysis of metallicities of Open Clusters is based on 
two different samples, which were constructed using differ- 
ent criteria to select the results from the literature. The first 
sample only includes [Fe/H] measurements obtained with 
high-resolution spectroscopy. This is basical l y the same sam- 
ple of 45 object collected bv lMagrini et all feoogh . updated 
with a few recent observations. These includ e NGC 6633, 
NGC 2539, NGC 24 47, IC 2714 and NGC 5822 dSantos et all 
l2009h . NGC 7160 iMonroe fc Pilachowskil (|201Ch , CR110. 
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Figure 2. The gravitational potential of the Galaxy with a cen- 
trigal barrier. Two curves are shown, corresponding to minimum 
energies (circular rotation) at 6 and 11 kpc. The radial position 
of circular rotation is indicated by a dot. The radial range of 
the stellar orbits for a small energy perturbation is indicated by 
dashed lines. 



NGC 2099, NGC 2024 and NGC 7 789 (jPancino et alJl20ld ). 
Tombaugh 2 dVillanova et alj|20ld l . Two measurements pre- 
sented by iGrattonl (|2000D . of IC4725 and NGC6087, al- 
though not so recent, were added to this sample. 

The second sample contains metallicity derived from 
spectroscopy as well as from photometry and other meth- 
ods. These data are less reliable, and their use is justified by 
the larger number of measurements available (131 compared 
to 58 for the first sample). This second sample is based on 
the on-line version of the New Catalogue of Opticall y visi- 
ble Op en Clusters and Candidates published by iDias et aU 
((2002J), and often updatecQ. This catalog contains a list of 
153 clusters with available metallicity measurement, and 
gives the adopted reference for each entry, as well as the 
method of observation and estimated error. The metallicity 
indicator used is [Fe/H], usually given in the original refer- 
ences. However, in a few cases the metallicity included in the 
catalog is Z, from which the Fe abundance ca n be recovered 
using the relation [Fe/H] = log (Z/0.019) (see[c arraro et al.l 
Il999h . 

Minor changes or updates were performed for a number 
of objects, often giving preference to metallicity and dis- 
tance from the same author. For NGC 2 168 w e adopted 
the distance given by IPohnl fe Paunzenl (l201fj|)(0 8 kpc) ; 
for NGC 2489, the distance given by iPiatti et all (|2009l ) 
(1.8 kpc); for NGC2335, the same distance used by TAAT 
(1.15 kpc); for NGC188, the distance given bv lMeibom et all 
1 20091) (1.77 kpc). and f or NGC2360, the distance given 
by lHamdani et alj i|2000l ) (1.09 k pc); for Berke l ey 20 and 
Collinder 261, we used the d ata of Sestito et al.l (|2008l ); for 
Berkeley 31 the distance by lYong et all l|2005l ) (5.3 kpc); 
for IC4756, NGC 2 682 and NGC 3680, the [Fe/H] data of 
ISantos et all l|2009l) . 

In the analysis of this second sample, we only make use 
of data with errors smaller than 0.2 dex in [Fe/H]. A num- 
ber of authors do not give an error estimate (see the on-line 
catalog), but we adopt an error estimate of 0.1 dex based 
on the sim i larity of the technique with other work. Data of 
lAnn et alj (|2002h are an exception, for which we believe that 
the error is larger than the limit adopted here, because some 
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discrepancies were found in two Op en Clusters that are in 
common with Hascgawa et all (|2004l ). and because for some 
clusters it seems that Z is given instead of [Fe/H], which is 
not clear. The data from Hasegawa et al., (not included in 
the on-line catalog), on their turn, are not included because 
the error estimated by the authors is 0.3 dex. The error on 
the metallicity of Berkeley 21 (the Open Cluster wit h the 
lowest metallicity) is about 0.3 dex i Tosi et al. 19981). and 
that of NGC6451 is about 0.5 dex i|Paunzen et al.ll2003l ). It 
should be noted that our procedure of eliminating objects 
with large errors on metallicity has affected more severely 
objects with [Fe/H] < -0.5; it may happen that our sample 
is biased in showing only few clusters with very low metal- 
licity. Note that there is a large overlap in the list of clusters 
between the high-quality sample (first sample) and this one. 

For both samples (high-resolution, and general) we re- 
computed the galactocentric distances of all the objects us- 
ing Ro = 7.5 kpc instead of the IAU recommended value 
8.5 kpc. This is a more realistic estimate of Ro, as for 
instance the (prefer r ed) M odels 1 and 2 from Table 2 of 
iMcMillan fe Binnevl (|2010l ) assume Ro = 7.3 and 7.8 kpc 
respectively. A higher value of Ro would imply an unre- 
alistically h i gh cir cular speed. See also the recent work of 
ISofue et alj (|201ll ). But the use of this shorter scale is not 
something new in the literature, since many specialists in the 
Galactic distance sc ale recognized tha t Rp is about 7. 5 kpc; 
see the reviews by iFich fe Tremaind (Il99ll) an d by iFeastJ 
(2008), or the papers by Glushkova et al.l ( 19991 ) and others 
of the same group. It is preferable to use this value of Ro 
when precise distances in different directions are to be com- 
pared, like for instance when stellar distances are compared 
with the position of resonances. 

In Figure 3 we present [Fe/H] as a function of galac- 
tocentric distance, for the sample of Open Clusters with 
high-quality measurements. The data of the second sample 
(lower quality but larger sample) are presented in a similar 
way in Figure 4. In both figures different colors were used to 
indicate young, intermediate age and old Clusters. In Figure 
5 the same data of Figure 4 (larger sample)is presented with 
histograms on the two axes. 

3.1 The step and the gap at corotation 

The first remark to be done is that the step in the Fe 
abundance first observed by TAAT is confirmed; it appears 
clearly in Figure 3 which presents high-quality data that 
were not available at the epoch of TAAT's paper. The height 
of the step is about 0.3 dex. The Galactic radius at which 
it is observed (8.5 kpc) is the same found by TAAT after 
correction for the different values of Ro- In addition to the 
step, one can see a gap in the distribution of clusters at 
Galactic radius 8.5 kpc, which is clearer in Figure 4 and in 
the lower histogram of Figure 5, due to the larger number of 
Clusters plotted in these figures. The radius of the gap coin- 
cides with corotation and with t he Cassini-like ring-s haped 
region void of gas discovered by lAmores et alj <|2009h . The 
gap in the distribution of clusters tells us that the star for- 
mation is hindered in the region of gas void, while the step 
in the metallicity distribution indicates that the gas void 
forms a barrier that has prevented the transportation of gas 
from one side to the other, so that the metallicity evolved 
independently on the two sides. 
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Figure 3. Fe abundance of Open Clusters normalized to solar as 
a function of radius measured by different authors. Only results of 
high-resolution spectroscopy are presented (our Sample 1). The 
errors on [Fe/H] are taken from the original references. Galato- 
centric distances were re-calculated for Ro= 7.5 kpc. Different 
colors correspond to ranges of ages of the Open Clusters; blue: 
age < 200 Myr, green: 200< age < 1200 Myr, red: age > 1200 
Myr. The dashed horizontal lines indicate the average metallicity 
on both sides of corotation. 




Figure 5. The same data of Figure 4, without error bars and 
without age information. On the two sides of the plot, the his- 
tograms of the data on the corresponding axis are presented. The 
green line indicates the corotation radius, which is a frontier be- 
tween two distinct regions of the disk and also corresponds to a 
minimum of star formation rate. The grey bands emphasize con- 
centrations of clusters at metallicities corresponding to peaks in 
the histogram on the left side. 
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Figure 4. Like Figure 3 for our Sample 2 of Open Clusters, which 
includes [Fe/H] results of different quality, based on spectroscopy, 
photometry and other methods. Measurements with errors larger 
than 0.2 dex were excluded. 



3.2 The flatness of the external gradient and 
cluster ages 

The metallicity distribution as a function of radius con- 
firms that the metallicity gradient is flat in the outer 
part of the disc, a fact w hich was observed by many au- 
thors, a mong them TA A T.ICarnev et all (120051). I Yo-ng et all 
(2005), ICarraro et all i|2007h and iMagrini et all (|2009). 



There was not yet a consensus on where the flat part be- 
gins, but from Figures 3 to 5 it seems clear that the radius 
at which the "base level" of metallicity of the outer disc 
starts is the corotation radius. 

In Figure 3 and 4 we separated the Open Clusters in 3 



groups of age, the young ones (age < 200 Myr), intermediate 
(200 < age < 1200 Myr) and old (age < 1200 Myr). One 
can see that there is no young cluster with [Fe/H] smaller 
than about -0.3 dex, on the inner side of corotation. The 
young clusters present are IC 2581, NGC 2343 e NGC 6716. 
Surprisingly, many metal-rich clusters with [Fe/H] about 0.1 
dex are observed in the same range of Galactic radius. This 
indicates that large metallicity differences are present in the 
gas situated at different azimuthal directions in a same in- 
terval of Galactic radius, which suggests that a fast mixing 
of the gas in the azimuthal direction, usually assumed in 
chemical evolution models, is a poor approximation. 

There is only one young cluster beyond corotation in 
Figure 4, Dolidze 25. This brings the question: at what 
Galactic radius does star formation stop? Star formation 
tracers like HII regions, compact HII regions revealed by CS 
emission and other tracers ar e seen up to about 14 kpc, e.g. 
see L+5 and IVazquez et all l|2008h . In Figure 3 and even 
in Figure 4 we are working with small samples of objects 
with well measured metallicity, so that we should not con- 
clude from these samples that young clusters do not exist at 
larger radii. 

In Figure 6 we examine the age versus Galactic radius 
relation of the larger sampl e of cluster s (126 6 objects ex- 
tracted from the web page of lDias et all (2002)) which have 
measured distance and age, without any other requirement 
(for this study the existence of high-resolution spectrometry 
is not so relevant). This figure confirms that young clusters, 
if they exist, are very rare beyond 14 kpc. The most distant 
object among the young clusters (in the lower part of the 
Figure 6) is Teutsch 45, at Galactic radius 14.3 kpc, which 
presents, however, a large error on distance (± 1.6 kpc). 
Other objects which were previously believed to have similar 
Galactic radii, like Dolidze 25 and NGC 1893, had their dis- 
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tance corrected based on rec ent publications l|Delgado et al.l 
l20ld ,l Prisinzano et alj|201ll ) . We note however that Carraro 
et al. (2010) report the detection of diffuse groups of young 
stars in the third quadrant of the Galaxy, at distances be- 
tween 14 and 19 kpc (converted to our short scale). The 
estimated distance of VdB-Hagen 4, which contains B stars , 
is of the order of 19 kpc. Similarly iMoitinho et all (|2006l ) 
detected " blue plumes" , which are groups of young stars up 
to distances of 16 kpc (in our scale) in the third quadrant of 
the Galaxy (the other "normal" clusters listed by the same 
authors are already included in the cluster database that we 
presented) . 

The flatness of the gradient of the outer disk (Figures 3 
and 4) deserves an explanation. One possibility is the flow of 
gas from corotation to the external regions, which is a conse- 
quence of the interaction of the gas with the spiral potential 
perturbation. The radial flow is predicted t heoretically and 
confirmed by hydr odynamical simulations (|Mishurovll2000l . 
iLepine et al.ll200lh . This flow transports the relatively high 
metallicity gas of the corotation zone towards outer regions, 
and flattens the gradients. 

Another possible reason for the flatness of the gradients 
is that a large fraction of the clusters that are presently ob- 
served at galactocentric distances of 15 to 20 kpc were born 
at smaller distances. This is naturally expected based on 
the discussion presented in Section 2. As we move towards 
larger radii, the curves representing the effective potential 
(Figure 2) become shallower, and for the same perturbation 
energy, the clusters cover a larger Galactic radius range. 
This hypothesis aims to explain why there are no (or almost 
no) young clusters beyond about 14 kpc. The conditions for 
star formation (probably, the gas density) are not fulfilled 
at large radii; the age of the clusters situated there is suf- 
ficiently large to allow them to be born in an inner radius 
and to have traveled to their present position. An additional 
argument in favour of the hypothesis that the outer clusters 
are travelers born in inner regions is their distribution in 
Galactic component z (direction perpendicular to the Galac- 
tic plane), shown in Figure 7. Simple integration of orbits 
of stars (or clusters) in the Galactic gravitational potential 
shows that if they are launched from a radius of 12 kpc, 
when they reach a radius about 18 kpc, their z coordinate 
can easily reach values of the order of 2 kpc, similar to the 
observed scale height. In such experiments the same initial 
velocity perturbations in the z direction required to explain 
the local (solar neighbourhood) scale height is used. Since 
the gravity of the disk decreases with Galactic radius, the 
amplitude of the oscillation in the z direction increases, if a 
star moves outwards. The good match of the scale-heights 
points towards a common origin of the velocity perturba- 
tions. The smooth increase of average age with distance seen 
in Figure 6 is also an argument in favour of the "travelling 
clusters" hypothesis. Note, however, that there are indica- 
tions that the two outermost old Open C lusters, Berkeley 29 
and S aurer 1, are of extragalactic origin l|Carraro fc Bensbvl 
I2009T ). 



3.3 A possible zone of enhanced star-formation 
efficiency between the 4:1 resonances 

Does the 4:1 resonance have any visible effect on the star- 
formation rate? This resonance, which is the ILR of the 4- 
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Figure 6. The age of Galactic clusters as a function of Galactic 
radius, for the larger sample of clusters which have estimated 
distance and age (with no need for metallicity measurements). 
The two vertical dashed lines are at symmetrical Galactic radii 
with respect to the Sun (at 6 kpc and 9 kpc) to emphasize that 
the cluster distribution is not symmetrically distributed around 
the Sun. 
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Figure 7. The distribution of Galactic clusters as a function 
of Galactic radius and of distance Z from the Galactic plane, 
for the larger sample of clusters which have estimated distance 
and age. The ages of the clusters are indicated by colors, like in 
Figures 3 and 4. The blue color symbols (young Clusters) are 
densely concentrated near the Galactic plane, but are hidden by 
the green symbols. Note that the distance scales are different on 
the two axes. The older clusters reach larger \Z\ and larger R 
values. 



arms structure, as previously discussed, is at 6.2 kpc, see 
L+5. Since our Galaxy has a strong 4-arms component, the 
frontier of the region where 4 arms are allowed to exist with 
the region with only 2 arms is in principle an important 
boundary. The star-formation machine (the spiral arms) and 
the mechanism that produces the inflow of gas have their 
power increased beyond that radius due to the larger number 
of arms. All the higher resonances, some of them possibly 
having associated spiral arms, are situated between the 4:1 
resonance and corotation. Indeed, both the lower histogram 
in Figure 5, and Figure 6, which exhibits a larger sample of 
clusters, show a strong increase in the density of clusters at 
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about 6.2 kpc. Of course, the interstellar extinction could 
be a cause of the cut-off in the counts of clusters at Galactic 
radius about 6 kpc. However, the 4:1 resonance is only 1.3 
kpc from the Sun, which corresponds to a visible extinction a 
little larger than 1 magnitude in the direction of the Galactic 
center. In the opposite direction, many clusters are observed 
at Galactic radius 12 kpc, which is much farther from the 
Sun (5.5 kpc). Therefore, the abrupt increase in the density 
of clusters beyond 6 kpc seems to be real, not an effect of 
extinction. It can be seen in Figure 6 that there are much 
more Clusters on the right side of the vertical line at 9 kpc 
than on the left side of the line at 6 kpc, although the two 
lines are equidistant from the Sun. The outer 4:1 resonance 
(fl p = fl+fc/4) is situated at 12 kpc. For the same reasons we 
could expect a decrease of the star formation rate and cluster 
counts beyond that radius. However, the frontier does not 
seem to be so well defined in this case. The probable reason 
is that the completeness of the sample of clusters is smaller, 
at distances of about 5.5 kpc. Nevertheless, it seems that 
the Galactic radius interval 6 to 12 kpc, the region between 
the two 4:1 resonances, except for the corotation gap, is a 
privileged zone for star formation. As we discussed, the old 
clusters situated at large distances possibly originated in the 
6-12 kpc zone, and also most of the intermediate-age clusters 
situated in the neighbouring region 12-14 kpc could possibly 
be born inside the 12 kpc radius. 

At 12 kpc there is some indication of a minor gap in 
the distribution of clusters (see eg the histogram in Figure 
5). In the alternative interpretation of the spiral structure 
discussed in the introduction, in which a secondary pattern 
speed coexists with the main one, 12 kpc is the radius of 
corotation of the outer pattern. Some of the properties of 
the distribution of clusters could also be explained by this 
interpretation. 

3.4 Fine structure of the metallicity distribution 

Another interesting aspect of the metallicity gradients 
shown in Figures 3 to 5 is that they present concentrations 
of clusters around a number of horizontal (constant metal- 
licity) lines. Some of the peaks of the histogram on the left 
side of Figure 5 seem to be significant. For instance the peak 
at 0.15 dex also corresponds to a group of Clusters with the 
same metallicity in the sample with only high-quality mea- 
surements shown in Figure 3. We prefer to postpone the dis- 
cussion on concentrations of metallicities around preferential 
values to the later section on Cepheids, since the metallic- 
ities of these stars, which form a completely independent 
sample, confirm our statement. We remark that in Figures 
4 and 5 several old clusters situated inside the corotation ra- 
dius have high metallicities ([Fe/H] « 0.15 dex), which are 
larger than the average value for young clusters. This poses 
the interesting question: is the metallicity in the solar neigh- 
bourhood decreasing with time (more recent clusters have 
smaller metallicities) or did the old clusters situated around 
the solar radius have born in inner regions of the Galaxy? 
Still another possibility would be to associate this observa- 
tion with the azimuthal differences in chemical abundances, 
that are disussed in a later section. This question shows 
again the importance of taking into account the radial mo- 
tion of the tracers to understand the metallicity gradient. 
In Figure 6 one can see a concentration of clusters with 



age 10 7 years; this can be interpreted as a rounding effect 
in the papers reporting ages, rather than an evidence for 
a star-formation episode. It should be remembered that the 
large sample presented in that figure (1266 objects) contains 
data with non-uniform quality. 

We will not discuss here the abundances of a-elements 
in Open Clusters, as only a small number of clusters have 
been measured. The a-elements can be better analyzed 
based on the sample of Cepheids, discussed in the next sec- 
tion. 



4 THE GRADIENTS BASED ON CEPHEIDS 

An impressive set of metallicity data for 130 Cepheids, 
with abundances derived for many elements based on 
high resolution spectra, is presented in the series of pa- 
by A ndrievsky a n d coll ab orators: lAndrievskv et ail 
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stitute the main source for the present analysis. The range 
of Galactocentric distances covered by the sample is 4 to 16 
kpc. 

We revised the distances of the Cepheids to the Sun 
given in those papers as next explained, and subsequently 
re-computed the galactocentric distances using Rq —7.5 kpc, 
as it was done for the Open Clusters. The stellar distances in 
the original papers on metalli city measuremen t s were usu- 
ally taken from the catalog of lBerdnikov et all (2000) with 
on-line updates. We re-estimated them using the expres- 
sion for intri n sic co lors as a function of period given by 
lAbrahamvanl l|2003h and in addition, using R„=Av/E(B- 
V)= 3.5 instead of the usual 3.1 used by Berdnikov et 
al.. The reason for this change is that R„ depends on the 
spectral type of the stars (a correction must be made to 
the effective wavelength of the B and V filters when they 
are convolved with the en ergy distribution of the stars (see 
lAzusienis fc Straizvslll969T ) . L+5 found that the distances of 
the Cepheids calculated this way were better correlated with 
the kinematic distances than the original distances given by 
Berdnikov et al. 

These minor corrections of distances do not change the 
main picture of the gradients that emerged form the series of 
papers of Andrievsky and collaborators, that we described 
in the Introduction. 

We present in Figure 8 the Fe abundance in Cepheids 
as a function of Galactic radius. In this plot we made small 
corrections to the values of [Fe/H] given in two papers. The 
papers of Andrievsky et al. and Luck et al. are based on the 
use of the same approach of spectroscopic analysis (the same 
system of oscillator strengths and the same methods of de- 
termination of the stellar fundamental parameters) so that 
the results are directly comparable. In spite of this, a mi- 
nor mismatch appeared in the average metallicity for a large 
number of stars in the same galactic radius (there were no 
st ars in common). W e added 0.03 dex to the measurements 
of lLuck et all (2006) to correct for the difference in the aver- 
age value of [Fe/H], considering the data of Andrievsky et al. 
as the correct ones. Anoth er correction , of lar ger magnitude, 
was applied to the data of lYong et alj (|2005l ). Their sample 
contains two stars in common with the series of Andreievsky 
et al: NT Pup and CE Pup, for which Yong et al. obtained 
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Figure 8. Fc abundance of Cepheids normalized to the solar val- 
ues as a function of Galactic radius measured by different authors: 
black dots, the seri es of papers b y Andriev sky and collabora - 
tors, orange dots .lLu ck et al l ||2006T), blu e dots lYong et all^OOSl) . 
and green dots, iKovtvukh et al' 1 720051) . Minor corrections have 
been applied to the data of some of these authors (see text). The 
horizontal dashed lines emphasize metallicities around which the 
Cepheids seem to be concentrated. The galatocentric distances 
are based on Rg = 7.5 kpc. 



[Fe/H] values lower by 0.25 and 0.2 dex respectively. The av- 
erage difference for the stars situated in the region of overlap 
of the two samples was 0.3 dex. We added 0.25 dex to the 
[Fe/H] results of Yong et al. to match the data of Andriesky 
et al. This correctio n has an uncerta i nty of the order of 0.05 
dex. The sample of IKovtvukh et all (|2005l ) contains 4 stars 
in common with Andrievsky et al.: RX Aur, RW Cam, T 
Mon, SV Mon; the average metallicity difference for these 
stars is negligible. 

At a first glance, the traditional inclined straight line 
used in the past to fit metallicity gradients as discussed pre- 
viously, would seem to be justified, to describe the ensemble 
of Cepheid data. However, a step function (or Heavyside 
function), not shown, with a step of 0.3 dex at 9.5 kpc, pro- 
duces a fit of the same quality (Chi-square divided by the 
number of degrees of freedom), if we restrict the fit to the 
interval 6 to 13 kpc. In other words, this plot cannot be used 
to dismiss the existence of a TAAT-like step in the sample of 
Cepheids. The magnitude of the step, of about 0.3 dex, seem 
to be a little smaller than the one found for Open Clusters, 
but the difference is within the errors. 

The gap in the distribution of Galactic radii of the 
Cepheids, similar to that of the Clusters, can be seen at 
about 8.7 kpc. It seems to be narrower and a little displaced 
with respected to the corotation radius (8.4 kpc). One ex- 
planation could be that the Cepheids are on the average 
older than the young Clusters and had more time to par- 
tially invade the corotation gap. The r ange of ages of th e 
Cepheids is about 50 to 300 Myr (e.g. iBono et al] 120051 ). 
Another possibility is that the errors in distance are larger 
for Cepheids. 

A number of horizontal lines are plotted in Figure 8 
at [Fe/H] values which seem to concentrate a large quan- 
tity of stars. Some of these lines coincide with the [Fe/H] 
values with relatively larger density of clusters in Figures 3 
to 5, like for instance at 0.0 dex and 0.15 dex. This is not 



surprising, since the two samples are constituted of popu- 
lations of young objects; both can be considered to repre- 
sent the chemical abundance of the interstellar gas at their 
birthplace at the present epoch, since 200 Myr (the order 
of magnitude of the ages) is a small interval of time for the 
chemical evolution of the disk. 

To produce the horizontal alignments of Cepheids seen 
in Figure 8, there must be relatively few regions in the Galac- 
tic plane with active star formation. In any one of these re- 
gions, stars are formed from the local gas with a same metal- 
licity, and then start their travel in radial distance from the 
Galactic center according to their initial perturbation veloc- 
ity. The gas of the disk possibly presents a smooth variation 
of metallicity with radius, except for the step at corotation. 
If there were star forming regions distributed all over the 
Galactic plane, with a smooth distribution of metallicity, it 
would not be possible to distinguish preferential metallici- 
ties in the distribution of clusters. Certainly the Cepheids 
are born in spiral arms, like the clusters (in the case of the 
clusters this was shown by direct i ntegration of their or - 
bits for a time equal to their age, see iDias fc Lepind i|2005h . 
Therefore, the spiral arms are the obvious candidates to be 
the discrete star-forming regions required to explain the ex- 
istence of peaks in the metallicity distribution. Then, a nat- 
ural question is: many arms are elongated structures, span- 
ning a range of Galactic radius, why don't they produce 
a broad distribution of metallicity that would smooth out 
the horizontal lines seen in Figure 8? One possibility is that 
some portions of the arms, like the corner of the square- 
shaped resonant arm discussed by L+5, are more efficient 
star formation regions than the others, and dominate. An- 
other possibility is that due to the flow of gas along the arms 
the gradients of metallicity are small within a given arm. 



4.1 The a-elements 

In order to reduce the errors on the abundance of a- 
elements, we averaged five elements. These are O + Si + 
S + Mg + Ca, taken with equal weights. We shall designate 
them as a-elements for simplicity. The choice of these ele- 
ments was based on the fact they are classical a-elements, 
and they were measured by different authors who performed 
spectroscopy of C epheids. Since O was not measured by 
lYong et al.l (|2005l ). for the data of these authors we calcu- 
lated the average of the remaining elements. The results are 
shown in Figure 8. We also performed minor corrections to 
the data presented in this figure, to match of the alpha abun- 
dance data of diffe r ent a uthors: the correc tions were -0.04 
dex fo r [Luck et all (|2006l ) and -0.1 dex for IKovtvukh et al.l 
i|2005l ). 

The step-like decrease in the abundance of a-elements 
at corotation looks similar to that of the Fe abundance. We 
plotted two horizontal lines in the Figure, to indicate the 
average metallicity inside corotation and the "base-level" 
beyond corotation. The difference between them is of 0.25 
dex, a little smaller than the step in Fe abundance but the 
difference is not significant compared to the scattering of the 
points. 



9 



1 .0 [ 



Cepheids 



0.5 - 



0.0 



15 - 



-1. 



8 10 12 14 16 

Galactic radius (kpc) 



1 .o r 



Cepheids 



0.5 - 



° 0.0 



-0.5 - 



8 10 12 14 16 

Galactic radius (kpc) 



Figure 9. Abundance of a-elements (average of O + Si + S + 
Mg + Ca) in Cepheids normalized to the solar value, in dex, as a 
function of Galactic radius (Rq= 7.5 kpc) measured by different 
authors. The colors are the same as in Figure 7. The dashed lines 
indicate the average abundance of a-elements in regions inside 
corotation (-0.03 dex) and the baselevel at large distances (-0.28 
dex). 

4.2 The [a/Fe] ratio 

The obtained [a/Fe] ratio, where a is composed of the same 
elements above, is presented in Figure 9. The same minor 
corrections that were applied independently to Fe and a- 
elements abundance were used to compute [a/Fe]. The ra- 
dial distribution shown in Figure 9 is surprisingly flat, show- 
ing an homogeneous abundance pattern all over the disk. 

Although our discussion here has its focus on Cepheids, 
we have to make comparisons with results for different ob- 
jects. Our_j«suh^se«m_to be in contradiction with those 
from ICarnev et al.l (|2005l ) , who have observed a relatively 
large (0.2 dex) increase of the [a/Fe] ratio beyond 10 kpc. 
The origin of this discrepancy is mostly due to the Fe abun- 
dance, found by Carney et al. to be abo ut -0.5 dex in the 
outer disk, while lAndrievskv et al.l (|02 d) derive about -0.2 
dex, not much different from TAAT (~ -0.3 dex). Supposing 
that our results based on Andrievsky's data are correct, the 
uniformity of the [a/Fe] ratio with Gala ctic radius i s an ar - 
gument against the hypothesis raised bv lYong et all l|2005l ). 
according to which the outer Open Clusters could be asso- 
ciated with mergers, and consequently, have very different 
forma tion histories from t hose born in our Galaxy. 

lAcharova et alj (|2010l ) presented in the introduction of 
their paper a detailed discussion of the problem posed by 
the similarity of the Fe and a abundances. According to the 
traditional view, most of the Iron is synthesized in type la 
supernovae (SNe), whose progenitors are stars with ages of 
several billion years. On the other hand, a-elements are syn- 
thesized in type II SNe, associated with massive short-lived 
stars. When star formation begins somewhere, the abun- 
dance of a-elements grows quickly, while the Fe abundance 
only starts to in c rease much later. However, more recently 
iMatteucci et al. | l|200fj ) developped a new formulation for 
the Type la SNe rate, in which up to 50 per cent of the 
total Type la SNe should be composed by binary systems 
with lifetimes as short as 10 8 yr. Taking into account that 
type II SNe also produce Iron, at least 2/3 of the Iron is 



Figure 10. [a/Fe] as a function of Galactic radius. The color 
scheme is the same of Figure 7. 

produced in a fast process. Consequently, the [a/Fe] ratio 
is not very sensitive to the time elapsed since the beginning 
of star formation, and is a good indicator of a recent star 
formation history. In this interpretation, the flatness of the 
distribution in Figure 9 is not so surprising. 

Still, it is possible to obtain some insight on the star 
formation history from this figure. One can see a number 
of stars with [a/Fe] of the order of 0.1 dex, in the range 
of galactocentric distances about 7 to 13 kpc. This is an 
indication that there has been intense star formation activity 
in recent times in this Galactic radius range. Interestingly, 
the stars with high metallicity seen around 6 kpc in Figure 
8 do not have a high [a/Fe] ratio. 



5 THE AZIMUTHAL GRADIENT 

We show in Figure 11 the azimuthal distribution of Fe abun- 
dance in the sample of Cepheids. This angle is the galacto- 
centric angle in the Galactic plane, defined to be 0° in the 
direction of the Sun and positive in the direction I = 90°. 
Obviously, if we could observe all the 360°, no average gradi- 
ent would be seen. In this example we restricted the range of 
Galactic radius to go from 7 to 11 kpc. Very similar plots can 
be obtained if we use the abundance of a-elements instead 
of [Fe/H], or if we change the range of radius. Note that if 
we transform the two extremities of the line segment fitted 
to the data in the figure into points on the Galactic plane, 
their distance is about 6 kpc, taking the average Galactic 
radius equal to 9 kpc. Therefore, the azimuthal gradient 
transformed to distance units would be of the order of 0.05 
dex/kpc, which is not negligible compared to the radial gra- 
dients found in the literature. For instance, iPedicelli et al.l 
(2009) have found an overall radial gradient of 0.05 dex/kpc. 

We attribute this local azimuthal gradient to the spiral 
structure. For example, there is an excess of on-going star 
formation in the second Galactic quadrant (£ — 90 — 180°) 
compared to the third one, due to the presence of the Perseus 
star forming region. This can be seen in maps o f the spiral 
struct ure based on young tracers, like those from lHou et al.l 
|2009l ) and L+5. 

The main purpose of presenting this azimuthal gradient 
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Figure 11. The Fe abundance of the sample of Cepheids as a 
function of the azimuthal angle. This angle is the galactocentric 
angle in the Galactic plane, defined to be 0° in the direction of 
the Sun and positive in the direction I = 90°. Only the Cepheids 
in the range of Galactic radius from 7 to 11 kpc are included. 



Figure 12. The Oxygen abun dance as a function of radius in 
M 83 (NGC5236) measured by l|Bresolin et alJl200Sh . The gala- 
tocentric distances are in units of R25, equal to 8.4 kpc according 
to the authors. The tracers are HII regions. 



is to call attention to the risk of simple interpretations of 
radial metallicity gradients, and to emphasize the role of non 
axis-symmetric spiral arms. Obviously the usual assumption 
in chemical evolution models of almost instantaneous mixing 
of the enriched gas in azimuth is a poor approximation. It 
must be kept in mind that within 1 kpc from corotation 
(which includes the solar vicinity), the differential velocity 
of the gas with respect to the spiral structure is so small 
that it takes more than 1 Gyr for the gas to complete a turn 
around the Galactic center, in the frame of reference of the 
spiral structure. 

5.1 The metallicity gradient step in M83 

Since there is a widespread concept that abundance mix- 
ing in the azimuthal direction (along a circle around the 
Galactic center) is very efficient, we would like to present 
an additional example in which the simple parameteriza- 
tion of abundance patterns as a function of Galactic ra- 
dius only is a poor one, and to show that there can be 
more than one "organized" (in contrast with simple noise 
or random distribution) metallicity at a given radius. Fig- 
ure 12 presents the gr adient of Oxyge n abundance in M83, 
reproduced from iBresolin et"aH 1.2009). Coincidentally, the 
step occurs at about the same distance than in our Galaxy. 
Where the step appears we can find in M83 "overlapping" 
gradients which resemble the ones we can see in Figure 3. 
We do not claim, however, that 8.5 kpc is the corotation 
radius in the case of M83, since there are many discrepant 
determinations of this corotation in the literature. 



6 A LOWER LIMIT TO THE AGE OF THE 
SPIRAL STRUCTURE 

We have shown that the step in chemical abundance first 
noted by TAAT and confirmed in the present work coin- 
cides with the corotation radius; the step is explained it 
by the independent chemical evolution on the two sides of 



the barrier constituted by the ring void of gas and by the 
"pumping out" mechanism. On each side of the barrier, there 
is exchange of gas between neighbouring Galactic radii due 
to winds of massive stars, explosions of supernovae, and gas 
transfer in the spiral arms, so that relatively flat gradients 
have established. The discovery of the nature of the abun- 
dance step provides a unique opportunity to contribute to 
the debate on the age of the present spiral pattern of the 
Galaxy, a parameter of major importance for the under- 
standing of the physics of the spiral arms. A concept of 
"transient spir al arms" has been supported by a numbe r 
of papers (e.g. ISellwood fc Carlberd Il98i . ISellwoodll201ll ). 
According to that hypothesis, which is based on numerical 
simulations, the lifetime of the arms is only of a few Galactic 
rotations. However, based on similar t echni ques, but with a 
larger number of particles, iFuiii et alj(|2010h . concludes that 
the lifetime of a spiral structure can be of up to 10 Gyr. 

The existence of the step in abundance means that the 
corotation resonance stayed at the same position for a long 
period. If the spiral structure were a transient phenomenon, 
changing its characteristics and its corotation radius at time 
intervals of the order of a few galactic rotations, the abun- 
dance step would not exist. We next discuss a lower limit to 
the age of the metallicity step. Supposing that a new spiral 
structure established itself at some time in the past, with its 
corotation radius at its present radius, then different metal- 
licities started to build up on both sides of the resonance. A 
minimum elapsed time can be estimated by supposing that 
the production rate of metals was negligible on the outer 
side (of course, if the production of metals were not small 
on that side, a larger time would be required to build up 
the same 0.3 dex step). Taking t he slope of Fe enrichm ent 
as a function of time presented bv lChiappini et all (|l997h in 
their Figure 1, applied to the inner side, in the extreme case 
in which the enrichment is unimportant in the outer side, to 
build a 0.3 dex difference requires at least 3 Gyr. 

This number, which is a lower limit, is in contradic- 
tion with the concept of transient spiral arms, in the shorter 
pattern life-time version of it. A rough estimation of 3 
Gyr as the time elapsed since the last major change in 
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the spiral structure of the Galaxy was already proposed by 
lAmores et all (|2009l l based on the time required to build up 
the slope of the metallicity gradient of the youngest plane- 
tary nebulae. 



7 CONCLUSIONS 

The abrupt step in the metallicity gradient at about 8.5 kpc 
(fo r Rp= 7.5 kpc , or 9.5 kpc for i?o=8.5 kpc), first reported 
bv lTwarog et all (|l997h based on analysis of Fe abundance 
in a sample of Clusters, is confirmed by the additional data 
collected since 1997, and in particular, by the high-resolution 
sp ectroscopy data. A s imilar conclusion was already reached 
bv lCarnev et al.l (|2005l ). The step is also present in the sam- 
ple of Cepheids that we analyzed, not only for Fe abundance, 
but for the a-elements as well. 

The position of this step coincides with the corotation 
radius of the main spiral structure and is easily explained 
by the presence of this resonance (the possibility of the ex- 
istence of a second pattern sp eed was also discussed). As 
shown bv lAmores et al] (2009) there is a ring-shaped void 
of gas at this ra dius, which was pre dicted by hydrodynamic 
simulations (e.g. ILepine et al 1 l200lh . This void is caused by 
the pumping out effect of the spiral structure that produces 
an inward flow of the gas in the inner regions of the Galaxy 
and an outward flow in the outer regions. Therefore, the 
Cassini-like gap is a barrier that isolates the gas of one side 
from the other, so that the chemical evolution of the two 
sides are independent. One must remember that chemical 
evolution models are models for the gas; stars are not usu- 
ally regarded as a mean of transporting metallicity from 
one Galactic radius to the other. The corotation barrier only 
acts on the gas, not on the stars; many stars have orbits with 
strong deviations from the circularity, and there is no imped- 
iment for them to cross the resonance; what is more, stars 
on circular or bits themselves are scattered across the coro - 
tation radius (|Sellwood fe Binnevll2002l ; ILepine et al.ll2003h . 
and thus naturally cross the chemical barrier. 

The counts of Open Clusters suggest that in the Galac- 
tic radius range 6.2 to 12 kpc, that is, between the ILR and 
the OLR of the 4-arms mode, but excluding a narrow gap at 
corotation, the star-formation rate is larger than outside this 
region. This is not an unexpected result, since our Galaxy 
has an important 4-arms mode, which is only permitted in 
this region. 

Close to corotation, on the inner side of it, one can 
find Open Clusters with differences as large as 0.4 dex in Fe 
abundance at a same Galactic radius, which is much larger 
than the errors of measurements. The large dispersion of 
Fe abundance in this region can be attributed to the fact 
that the stars have been originated from different Galac- 
tic radii or different star-formation regions, through either 
non-circularity ("blurring") or angular momentum exchange 
of circular orbits ("churning"), making the variety of birth- 
places a more important cause of dispersion of metallicity 
than variations in age, except for the very old clusters. 

A fine structure study of the metallicity distribution in 
populations of clusters and Cepheids situated in a Galac- 
tic radius range of width 1-2 kpc, in any of the two sides 
of corotation, shows that the pattern of this distribution 
is not smooth but presents a number of peaks, shown by 



horizontal lines in Figure 5. Our interpretation is that the 
star formation process is active in a limited number of re- 
gions of the Galactic disk; each active region has its own 
well defined metallicity, which is shared by the stars born 
there. Due to their initial velocities at the instant of birth, 
the stars (or clusters) have non-circular orbits and usually 
present variations in their Galactic radius of the order of 2 
kpc in periods of time of the order of 100 Myr. In the ex- 
ternal parts of the Galaxy the amplitudes of the radial vari- 
ations get larger, enhancing the "blurring" effect, but the 
time-scales are also larger. The simultaneous presence, at a 
given Galactic radius, of stars that originated from different 
star-forming regions can be explained by the overlapping of 
orbits. The concept is illustrated by L+5, where we see that 
orbits of different natures can cross each other. In the same 
paper, the complexity of the local spiral structure is shown, 
with the young stars of the solar neighborhood distributed 
in spiral arms with many different orientations. 

In other words, the recent chemical enrichment of the 
Galactic disk seems to have its origin in a clumpy distribu- 
tion of dominant star-formation centers. As the stars move 
out from the star-formation centers, and present variations 
of their Galactic radius, they form the horizontal concen- 
trations of stars with similar metallicity like those seen in 
Figure 5. Certainly the star formation centers belong to the 
spiral arms. 

To further convince the reader that one cannot neglect 
the spiral structure in the interpretation of the chemical gra- 
dients, we showed that there is a significant azimuthal metal- 
licity gradient in a broad Galactic radius range around the 
Sun. With the same purpose, we showed an extragalactic ex- 
ample of the existence of a step in the gradient, associated 
with overlapping gradients (two different metallicity curves 
at a same radius). Analyses of metallicity gradients in our 
Galaxy or in external galaxies that are restricted to the de- 
termination of a single slope as a function of radius are poor 
ones. 

As a conclusion, we thus stress that several difficul- 
ties exist in the comparison of chemical evolution models 
with observations. Only recently, such models predicting 
the metallicity distribution of both th e stars and gas in a 
realist ic manner have been developed Schonric h fe Binnevl 
(2009), by taking into account both the angular momen- 
tum exchange across corotation of minor transient spiral fea- 
tures and the range of radii swept by stars on non-circular 
orbits. However, such recent models are still axisymmetric 
in essence, and do not take into account some of the basic 
features of our non-axisymmetric Galaxy, like the chemical 
barrier at the corotation of the main grand-design spiral pat- 
tern, the flow of gas with inverse direction on both sides of 
it, and a more correct expression for the star-formation rate 
taking into account the rate at which the spiral arms (where 
stars are born) are fed with gas. Consequently they cannot 
at present reproduce features such as the metallicity step 
and the overlapping abundance gradients observed around 
the corotation of the main spiral pattern, or such as the az- 
imuthal metallicity gradient observed in the vicinity of the 
Sun. More elaborated chemo-dynamical models taking into 
account the main effects of the grand-design spiral pattern 
of the Milky Way are thus now required to reproduce such 
observations. 
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